% 导入题目数据
data = [311   460   975   463   708   666   398   771   532   474
   538   740   651   458   407   420   467   207   457   337
   759   488   509   486   539   218   715   509   647   565
   314   613   530   578   599   319   574   647   730   481
   597   589   628   132   316   601   484   440   372   477
   497   591   243   587   172   668   865   362   678   382
   389   673   749   836   468   384   548   643   563   526
   749   487   417   649   570   214   527   308   553   743
   747   619   656   525   372   607   620   726   379   605
   280   586   763   851   653   492   528   607   590   590
   779   576   651   249   560   723   927   449   644   325
   619   734   320   599   754   433   521   971   175   582
   549   549   375   802   256   557   529   678   567   656
   627   502   708   531   503   452   677   524   539   212
   309   573   673   398   408   592   447   463   415   594];
data = data(:); % 转换成为向量

parm = wblfit(data); % Weibull 分布参数估计 parm(1) 为参数a parm(2)为参数b

p = 1/mean(data) * (1/9); % 非刀具故障的概率

F1 = @(x) wblcdf(x,parm(1),parm(2)); % 刀具故障分布函数
F2 = @(x) 1 - (1 - p).^x; % 非刀具故障分布函数

F = @(x) 1 - (1 - F1(x)) .* (1 - F2(x)); % 故障分布函数
P = @(i) F(i) - F(i - 1); 

tic
aveCost = objFun(18,316,100000,parm,p);
toc